Nonlinear Landau- Zener Tunnelling in Coupled Waveguide Arrays 
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The possibility of direct observation of Nonlinear Landau- Zener tunnelling effect with a device con- 
sisting of two waveguide arrays connected with a tilted reduced refractive index barrier is discussed. 
Numerical simulations on this realistic setup are interpreted via simplified double well system and 
different asymmetric tunnelling scenarios were predicted just varying injected beam intensity. 
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Landau- Zener tunnelling effect [T] has been first pro- 
posed for interpretation of atomic level mixing in predis- 
sociasion process [2]. More recently this model has been 
applied to explain transitions between Bloch modes in 
periodic systems, particularly for Bose-Einstein Conden- 
sates (EEC) [3, 4] and acoustic waves in layered and elas- 
tic structures [5] . Later on the same effect of Bloch mode 
transitions has been investigated in optical systems with 
variety of architectures: waveguide arrays with a step in 
a refractive index [6^, arrays with applied temperature 
gradient [7], curved waveguides [8], nematic crystals [9| 
and two dimensional photonic lattices [10], among others. 

The nonlinear extension of Landau-Zener model is first 
analyzed theoretically [III, [lH in case of EEC in op- 
tical lattices and asymmetric transition processes have 
been found in contrast to the linear limit. Later on 
this nontrivial behavior has been experimentally con- 
firmed [13] and generalized Landau-Zener transition for- 
mula has been analytically derived [is]. In all of these 
previous studies the tunnelling processes between differ- 
ent Bloch modes are considered, while, in principle, the 
same Landau-Zener tunnelling effect should take place 
for double well system. Indeed, very recently the tun- 
nelling scenarios between two spatial modes of EEC has 
been proposed [15|, JJi]. In the present paper we pro- 
pose to use two coupled waveguide arrays with a tilted 
(with respect to the waveguide direction) reduced refrac- 
tive index barrier for visual observation of asymmetric 
nonlinear tunnelling effects between the arrays. 

The refraction index profile of a suggested experimen- 
tal device is presented in insets of Fig. 1 and refractive 
index pattern could be realized either by microfabrica- 
tion [13 or by laser beams in photonic lattices p^]. The 
beam is injected either into the left or right array and 
its intensity has a harmonic profile across the injection 
array (this is schematically represented as black arrows 
with different lengths in the insets). In case of small in- 
tensity (linear regime) the light injected into the left ar- 
ray tunnels to the right array and vice versa. Increasing 
the injected beam intensity the symmetry breaks down, 
particularly, injecting the beam at the left it again tun- 
nels to the right array, while injecting the beam with the 
same intensity at the right it stays trapped there (see the 
corresponding graphs in Fig. 1). Further increase of the 
injected intensity leads to the beam trapping irrespective 
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FIG. 1: (Color online) Graphs (a) and (b): results of numeri- 
cal simulations (stationary spatial light intensity distribution) 
on the initial model ([Tj) with the amplitude 0.0018 of injected 
light wave pattern ^(F, 0). Insets show the device schematics 
and corresponding injection positions. 



to the place of the injection of the beam. Below I shall in- 
terpret this nontrivial effect via Nonlinear Landau-Zener 
tunnelling in a simple two degree of freedom system. 

Let us start from writing a wave equation for linearly 
polarized electric field in paraxial approximation: 



A(r,Z)^ + 1^1^^ = 0, 



(1) 



where I assume light propagation in nonmagnetic optical 
film (YZ plane) along Z direction, ^(y, Z) stands for a 
complex wave envelope, A(F, Z) = {n — no)/no is a linear 
refractive index variation along Y and Z axis, the dimen- 
sionless spatial variables Y and Z are scaled in units of 
inverse carrier wavenumber l//c, and this wavenumber is 
defined as /c = uquj/c with 00 being laser beam frequency 
and focusing Kerr nonlinearity is scaled to unity. Let us 
consider pinned boundary conditions (i.e. ^(F, Z) = 
for Y = Ymin-,ymax)i whilc periodic modulations to- 
gether with tilted reduced refractive index barrier (see 
insets of Fig. 1) are modelled as follows: 

A(r, Z) = 6ni sin^ {KY) - (5n2sech[(Y - rZ)/A], (2) 



2 




FIG. 2: (Color online) Reduction of the initial problem ([T]) 
with the refractive index profiles presented in the insets of Fig. 
1 to the effective double well system with "moving" barrier. 
The curved lines describe the light intensity dynamics for the 
linear (left graph) and asymmetric nonlinear (right graph) 
regimes. 



6ni and Sn2 stand for the amplitudes of periodic (scaled 
by K) modulation and barrier, respectively, while F and 
A define a tilt angle and a width of the reduced refractive 
index barrier. In numerical simulations on the model ([1]) 
the following values are fixed: 6ni = 0.002 and Sn2 = 
0.003. 

As well established, the problem of periodic array of 
effective waveguides could be simplified via tight binding 
discretization procedure [HI, [HI when one can present 
the envelope wavefunction "^{y^z) as an expansion over 
approximate gaussian eigenmodes of individual waveg- 
uides: 

^(r,Z) =^£j{Z)(l)j{Y) (l)j{Y) - e-^V^i'^-Ro)" 



where j numbers waveguide center positions and thus 
Rj = 7r(2j + 1)/2K. Then one gets a Discrete Nonlin- 
ear Schrodinger (DNLS) equation representation of the 
problem: 



■j+1 



,)-Vij,z)£,+x\£j\%=0, (3) 



where coupling constant C is calculated from the over- 
lapping integrals between neighboring waveguide eigen- 
modes, while X counts only the nonlinear overlap inte- 
gral of the single eigenmode. Effective potential barrier 
V(j, Z) could be approximated as (see Ref. (2Qj): 



V{j, Z) = Fosech [(TTj - rKZ)/KA] , 



(4) 



where Vq ~ 5n2 is a potential barrier height. Further 
reduction of (|3|) is made associating j with a continuous 
new spatial variable y = YK/tt and then defining z = 
CZ and rescaling £ £\fxlC exp[zCZ], V V/C one 
gets Nonlinear Scrodinger (NLS) equation in an external 
double well potential with "moving" barrier: 



.d8_ 

dz 
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F(^,z)f = 0. 



(5) 



The meaning of discretization of initial equation ([T]) 
getting DNLS ([3]) and then subsequent continuous ap- 
proximation to NLS equation (|5]) is that we get rid of pe- 
riodical modulation of refractive index which are present 
in initial equation ([T]). Indeed, in NLS (j5j) one is left only 
with "moving" barrier potential and the problem is re- 
duced to the two coupled waveguides case (see Fig. 2). 
Our aim is to give analytical consideration of the latter 
problem and present the interpretation of the numerical 
simulations undertaken on the initial equation ([1]). 

For the clarity of presentation let us choose symmet- 
ric boundaries —L < y < L and —io < z < i and we 
shall find the solutions of (|5|) in case of pinned boundary 
conditions S{—L,z) = £{L^z) = 0. Moreover, we will 
require small value for the barrier tilt parameter F ^ 
and then expanding expression for the barrier (j4]) over 
small parameter Tz we get 



V{y,z) = Vs{y)^TzVt{y) 



(6) 



Thus the potential barrier expression is split into sym- 
metric Vs{y) and antisymmetric Vt{y) parts with respect 
to the inversion transformation y ^ —y: 



Vs{y) 



Vo 



cosh[7rt//JsrA] 



Vosinh[7rt//i^:A] 



and let us build the stationary solutions of (j5j) with po- 
tential Vs (y) considering the second term Vt (y) as a per- 
turbation. Particularly, in the zero approximation we are 
left with the problem of symmetric double well potential 
which has two lowest eigenvalue symmetric ^~^{y) and 
antisymmetric ^~{y) orthonormalized solutions and one 
can simply construct from them two functions 



($+ + $-)/V2, 



($+ _ $-) /v/2 (8) 



localized at the left and right wells, respectively. 

Then one can separate the variables in £{y,z) estab- 
lishing a dimer model as 



£{y,z) = tl;i{z)^i{y) ^ tlj2{z)(l)2{y), i^l 



Pt (9) 



and substituting this into (|5]) where potential function is 
taken in the form (|6]), multiplying on ^i(^) and ^2(2/), 
then integrating over y and discarding common phase 
variables we recover the nonlinear Landau- Zener model 
PH in its standard form: 

= -azilj2 ^ viIji -r - |V^2p) V^2. (10) 

where the parameters could be calculated as follows: 

L L 

2r = I dyctA v= / (i#i,2 f ^ttt^ ~ '^s(y)(/>2,i ) 



dy'^ 
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FIG. 3: (Color online) Graphs (a) and (b) display dependen- 
cies of the propagation constant /3 and beam intensity at the 
left array iV^il^ on the parameter az according to equation 
(|lip in linear regime. Solid circles indicate initial states of 
the system and open circles stand for the final states in adi- 
abatic regime. Lower graphs display the results of numerical 
simulations on the initial model equation ([T)) when small in- 
tensity light is injected into the left (graph c) and right (graph 
d) arrays, respectively. 



and OL = — r (i^^i, 2^(^)01, 2 is an effective "acceler- 
ation" parameter. 

In the linear limit, i.e. when total intensity ipf -\-ip2 — 
Pt ^ the above equation is just an ordinary Landau- 
Zener tunnelling [1] which has a well known result that 
if "acceleration" a is large or coupling v is small the 
light remains in the array where it was injected initially. 
Otherwise, in adiabatic limit, the light tunnels to other 
waveguide array, thus the picture is symmetric. In our 



numerical simulations on the initial model equation ([T]) 
the barrier tilt angle F ~ a is small guaranteing adia- 
baticity of the process and one should see tunnelling of 
the light irrespective to the beam injection place. 

In order to clarify the behavior of the nonlinear two 
degree of freedom system (p!Q|) let us seek for the station- 
ary solutions in the form ?/^i^2 = |'0i,2| exp(z/3z) getting 
thus a quartic equation for the propagation constant /3 
and dependencies of intensities IV^ip and |?/;2p versus az: 



2|Vi 



rPt 1 



az 



(11) 

|^2|2 = rPt-|Vi|^ 



From these it automatically follows that the equation for 
j3 has four real roots if effective nonlinearity exceeds a 
coupling strength rPt > otherwise it has two real roots 
for fixed z (for more details see Ref. [ijj, ll2j). The cor- 
responding dependencies of /3 and on z for various 
total intensities are displayed in Figs. 3, 4 and 5. 
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FIG. 4: (Color online) The same as in Fig. 3 but with differ- 
ent parameters: Upper graphs correspond to the intermediate 
range of effective nonlinearity rPt = 1.5. Lower graphs dis- 
play the light intensity distribution in numerical simulations 
on initial equation ([T)) with wave envelope amplitude 0.0016. 



FIG. 5: (Color online) Upper graphs correspond to the strong 
effective nonlinearity range rPt > 2 and lower graphs display 
the light intensity distribution in numerical simulations with 
corresponding value of wave envelope amplitude 0.002. 

Let us first consider the linear regime rPt (see 
Fig. 3). Then the propagation constant is characterized 
by two solutions for fixed z. The solid circles in upper 
panels of Fig. 3 correspond to the starting point z = —i 
and we choose a£ = 2, when almost all the intensity is 
concentrated in the left waveguide array. This means 
that we should follow blue curve and in case of adiabatic 
process the system passes the point z = and propa- 
gation constant follows further the blue curve and as a 
result the system will end up with almost zero intensity 
at the left waveguide array, i.e. all the intensity should 
go into the right one. The same happens if the light is 
injected at the right waveguide array (in this case the 
system is initially on the black curve) and the light tun- 
nels into the left array. Indeed, in numerical simulations 
(see bottom panel of Fig. 3) on the initial equation ^ 
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FIG. 6: Numerical simulations on the discretized model Q 
with the parameters C — lOOmm"^ and x — 6.5mm~^VK~^. 
The and height of the barrier and its width are Vmax — 
0.6mm~^ and A = 25/im. Dashed lines indicate position of 
tilted reduced refractive index barrier and peak intensity in 
both cases is \Sj{Z = 0)\'^ax — 31W. 

it appears that beam injection with harmonic profile and 
small amplitude into the left array leads to the tunnelling 
to the right array and vice versa. 

Next let us consider the case when effective nonlin- 
earity exceeds a coupling constant rPt > v. First of all 
we note that as seen from the equations (pTj) a simple 
rescaling of all quantities with respect to the coupling 
constant v is possible and thus in our further analysis we 
can set V = 1 without loss of generality. Thus in case 
of effective nonlinearities rPt > 1 the propagation con- 
stant versus az diagram acquires a butterfly structure 
(see upper panels of Fig. 4) and starting again from the 
point az = —2 and with almost whole intensity in the 
left waveguide array (IV^ip — rPt) the evolution along z 
follows the blue line (see both graphs of the upper panel 
of Fig. 4) then it passes to the red line region where the 
butterfly structure begins [we note that different colors 
in Fig. 4 corresponds to the four different solutions of 
quartic equation (pT]) ]. By end of the red line there ex- 
ists a discontinuity, thus the system has to jump either to 
the blue line or to the black one even in ideally adiabatic 
case. On the other hand, starting with zero intensity 
at the left waveguide array (all the intensity is concen- 
trated at the right, i.e. IV^ip — 0) the system follows 
the black line and continues safely until the end point 
az = 2. This is the reason of asymmetric behavior of 



the nonlinear case, thus injecting the beam into different 
waveguide arrays one can end up at the same array. Dif- 
ferent amplitudes of initial harmonic beam profiles have 
been checked for this tunnelling scenario. For instance, 
in numerical simulations displayed in lower panels of Fig. 
4 the beam amplitude is 0.0016, while in Fig. 1 the am- 
plitude is 0.0018, and the asymmetric tunnelling regimes 
are observed. 

Increasing further the effective nonlinearity (in graphs 
of Fig. 5 we have chosen rPt = 2) the butterfly structure 
enlarges and the symmetric behavior restores. Particu- 
larly, and this is clearly seen from the both graphs of 
the bottom panel, injecting the light in the left or right 
waveguide array it remain trapped there. 

In order to make some predictions for realistic waveg- 
uide arrays the most convenient way is to consider dis- 
cretized equation ([3|). Then one is able to unify all the 
variety of refractive index profiles via a single coupling 
parameter C. In the numerical simulations presented in 
Fig. 6 I choose a coupling constant C/2 = Amm~^ and 
nonlinearity parameter is taken x = 6.5mm~^W~^ like 
in Ref. [17]. Assuming lattice spacing equal to 6/im the 
tilt angle of reduced refractive index barrier is chosen 
0.25mrad and then the samples of the length 60mm will 
be sufficient to see the effect of asymmetric tunnelling. 
As it was mentioned above the main reason of appear- 
ing of such asymmetry is a nonlinearity induced butter- 
fly structure in the reduced two degree of freedom model. 
In the realistic numerical simulations on discretized equa- 
tion (|3|) the peak intensity threshold for appearing of such 
a structure is \Ej{Z = 0)|^^^ = 291^ which is experimen- 
tally easily accessible [17] quantity. Above this intensity 
threshold the asymmetric tunnelling behavior takes place 
even for fully adiabatic processes, i.e. for very small tilt 
angles of reduced refractive index barrier, but then long 
waveguide arrays will be required in order to see the ef- 
fect. 

Concluding it could be stated that visual observation 
of nonlinear extension of Landau-Zener tunnelling in op- 
tical waveguide arrays has been proposed. All analytical 
predictions followed from simple two degree of freedom 
system are completely confirmed by numerical simula- 
tions on the model equations ([T]) and (|3]). 
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